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ABSTRA.CT 

A  study  of  the  parameters  involved  in  electron 
scattering  experiments  with  gaseous  targets  at  the  Naval  Post- 
graduate School  Linear  Accelerator  v/as  conducted.   The  para- 
meters included  gas  densities,  effective  target  thicknesses, 
and  radiative  corrections  to  experimental  cross-section  cal- 
culations .   To  compare  the  NPGLINAC  with  other  accelerators , 
an  electron  scattering  experiment  was  performed,  using  a 
hydrogen-helium  gas  mixture  target,  at  values  of  q^  from 
0.10  to  0.55  F~2.   The  resulting  charge  form  factor  for  he- 
lium was  fairly  consistent  with  previous  experimental  models.; 
however,  it  was  concluded  that  experiments  using  separate 
gas  targets  yield  greater  accuracy  than  those  using  a  gas 
mixture . 
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I.  INTRODUCTION 

A.   EXPERIMENTAL  OBJECTIVES 

The  electron  linear  accelerator  at  the  Naval  Postgraduate 
School  (NPGLINAC)  has  been  in  use  since  1966,  and  has  been 
described  in  theses  by  Barnett  and  Cunneen  [Ref.lJ  and 
Midgarden  [Ref.2].   Energy  loss  and  nuclear  structure  exper- 
iments have  been  conducted  using  solid  targets.   In  19  70,  the 
scope  of  experimentation  was  broadened  to  include  gas  targets 
in  addition  to  solid  targets.   While  the  technique  of  elec- 
tron scattering  measurements  remained  the  same,  there  V7ere 
several  nev/  problems  to  be  considered  v/hen  using  gas  targets 
instead  of  solid  targets.   Some  of  these  problems  included 
the  following  mechanical  considerations : 

(a)  Containers  had  to  be  designed  to  hold  gases  under 
high  pressures. 

(b)  A  method  of  cooling  the  target  area  and  keeping  it 
at  a  constant  temperature  had  to  be  designed  and 
installed. 

(c)  A  method  of  feeding  and  evacuating  the  gas  target 
chambers  had  to  be  devised,  along  with  a  method  to 
measure  the  pressure  inside  the  chamber. 

The  problem  associated  with  high  pressures  and  low  tempera- 
tures arose  because  it  was  necessary  to  maintain  the  gases 
at  high  densities  in  order  to  obtain  reasonable  counting 
rates  during  experiments.   Gases  at  atmospheric  pressure  and 


temperature  produced  counting  rates  that  were  very  low  and 
much  more  time  would  have  been  required  for  a  run. 

These  mechanical  problems  were  solved  by  Professors 
Bumiller  and  Buskirk,  Captain  Louis  Gaby,  and  the  NPGLIN/^C 
technicians.   A  diagram  depicting  one  of  the  gas  target 
chamber  arrangements  is  given  in  Figure  1. 

Apart  from  the  mechanical  problems  encountered  when  gas 
targets  were  introduced,  several  experimental  problems  also 
were  evident.  The  experimental  considerations  included  the 
following: 

(a)  Effects  of  the  target  chfimber  windows  on  electron 
scattering  experiments  had  to  be  determined. 

(b)  The  effective  target  thicknesses  of  gases  had  to  be 
determined;  i.e.,  only  electrons  scattered  from  a 
small  part  of  the  total  volume  of  gas  in  the  target 
chamber  made  their  way  into  the  spectrometer 
entrance. 

(c)  The  densities  of  the  gases  at  different  pressures 
had  to  be  determined  for  use  in  cross  section 
calculations . 

These  and  other  experimental  considerations  have  been  applied 
to  the  gas  target  experiments  conducted  at  the  Naval  Post- 
graduate School.   So  far  pure  deuterium  and  hydrogen-deuter- 
ium mixtures  have  been  studied  by  Gaby  [Ref.3]  and  Mader 
[Ref.4].   Hydrogen-helium  mixtures  have  also  been  studied  and 
are  reported  in  this  thesis. 
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The  purpose  of  this  paper  is  to  provide  guidelines  for 
use  in  further  studies  of  gas  targets  at  NPGLINAC.   Data 
from  deuterium  and  helium  experiments  were  used  in  preparation 
of  this  thesis. 

B.   EXPERIMENTAL  C7VLCULATI0NS 

One  of  the  most  important  considerations  in  a  cross- 
section  experiment  such  as  a  gas  target  study  is  the  evalua- 
tion of  radiative  corrections.   They  are  of  three  types:  the 
ionization  correction,  the  Bethe-Heitler  correction,  and  the 
Schwinger  correction.   These  corrections  are  described  in 
detail  by  Gordon  [Ref.5].   The  experimental  problems  of  tar- 
get chamber  window  thicknesses,  effective  gas  target  thick- 
nesses, and  gas  densities  all  appear  in  the  Bethe-Heitler 
correction.   It  v/as  therefore  considered  constructive  to  study 
these  parameters  and  report  their  significance  on  experimen- 
tal outcomes.   Likewise,  the  Schwinger  correction  was  studied 
for  similar  reasons.   The  Schwinger  correction  has  many 
forms  and  it  was  necessary  to  decide  which  form  was  to  be 
used  by  students.   Conformity  is  important  because  typical 
Schwinger  corrections  are  of  the  order  of  10-20%. 

In  addition  to  the  study  of  the  effects  of  experimental 
parameters  on  radiative  corrections,  a  cross-section  exper- 
iment was  performed  using  a  hydrogen-helium  mixture  as  a 
target.   From  the  relative  cross  sections,  a  calculation  of 
the  helium  charge  foim  factor  as  a  function  of  the 
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four-momentum  transfer  squared  was  performed.   The  results 
were  compared  to  those  obtained  at  Stanford  in  19  6  7  iRef.GJ 
and  to  the  Gaussian  model  of  the  form  factor. 

II.  RADIATIVE  CORRECTIONS 

A.   BETHE-HEITLER  CORRECTION 

The  Bethe-Heitler  radiative  correction  v/as  developed  to 
account  for  part  of  the  tail  below  a  peak  in  an  electron 
scattering  experiment.   After  incident  electrons  have  scat- 
tered off  a  target  nucleus,  they  lose  energy  in  the  form  of 
bremsstrahlung  as  they  pass  other  nuclei  and  the  atomic  elec- 
trons of  the  target  nucleus  itself.   The  less  energetic 
electrons  are  counted  at  much  lower  energies  than  the  peak 
energy,  theoretically,  all  the  way  to  zero  energy.   It  is 
not  possible  to  set  the  spectrometer  all  the  way  to  zero 
energy  because  of  the  length  of  time  involved  (the  ten- 
channel  system  can  only  monitor  an  energy  range  of  about 
three  Mev  per  setting)  and  the  background  radiation  levels. 
The  Bethe-Heitler  correction  is  an  estimate  of  the  number 
of  electrons  that  should  be  in  the  radiative  tail  due  to 
bremsstrahlung. 

Of  the  different  forms  of  the  Bethe-Heitler  correction 
term,  the  one  used  was  developed  by  Tsai  [Ref.7]  and  appears 
in  a  thesis  by  Gordon  [Ref.5]. 
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^B  =   I  K  Tiw  +  1/2  bt]  In  (E^^/n^AE)  + 

K  ^fw  +  1/2  bt]  In  (E3/AE)]  (2-1) 

where  the  b's  are  functions  of  the  atomic  number  g  only, 
T-L^  =  entrance  windov7  thickness 
Tf^  =  final  window  thickness 

t  =  effective  target  thickness 
Ej^  =  incident  electron  energy 
E3  =  elastic  peak  energy 
AE  =  E3  -  (lower  limit  of  electron  spectrum) 

and  n    =•  recoil  factor  =  1  +  2E2^  SIN^  6 

Mc2       2 
where  M  =  mass  of  target  nucleus 

and    6  =  scattering  angle. 

The  Bethe-Heitler  correction  is  applied  to  the  differen- 
tial cross  section  in  the  form  Kg  =  e~  B.   in  a  ratio  experi- 
ment using  a  two-gas  mixture,  a  relative  cross  section  is 
determined,  and  the  ratio  of  Kg  for  the  two  gases  is  enter- 
ed.  The  ratio  is  very  nearly  equal  to  one  in  most  cases. 
The  Bethe-Heitler  term  becomes  much  m^ore  significant  in  a 
pure  gas  experiment.   For  the  hydrogen-helium  runs  the  ratio 
of  the  corrections  was  taken  as  unity. 

In  the  Bethe-Heitler  correction  term  [Equation  2-1]  the 
important  experimental  parameters  include  the  window  thick- 
nesses and  the  effective  tarcet  thicknesses.   For  this 
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experiment  the  window  thicknesses  of  the  target  chambers 
were  knov/n  (1.0  mil  stainless  steel  for  the  two-inch  dia- 
meter chamber  as  well  as  for  the  three-inch  diameter  cham- 
ber) .   The  effective  target  thicknesses  were  determined 
from  the  experimental  cross  sections,  using  the  assumption 
that  the  counting  system  was  100%  effective.   The  calcula- 
tions involved  the  density  of  the  gas  in  the  target  chamber. 
The  gas  density  also  appears  in  the  ionization  correction 
term.   The  virial  equation  of  state  was  used  to  calculate 
gas  densities,  to  take  into  account  deviations  from  the  ideal 
gas  law.   The  virial  equation  of  state  and  the  method  of 
calculating  effective  target  thicknesses  are  described  in 
Appendixes  B  and  C. 

The  effective  target  thicknesses  were  caicuiatea  for 
hydrogen  and  deuterium,  using  data  from  the  helium  runs  and 
Mader's  deuterium  runs  [Ref.4].   The  results  are  shovm  in 
Table  I. 

It  was  noted  from  the  results  that  effective  target 
thicknesses  were  somev/hat  dependent  on  energy  and  very 
dependent  on  angle.   The  effective  target  thickness  v;as  a 
minimum  at  6  =  90°,  and  increased  in  proportion  to  the  devi- 
ation of  the  angle  from  90°.   For  a  particular  run,  the 
effective  target  thicknesses  of  the  two  gases  should  have 
been  equal.   Discrepancies  of  from  less  than  1%  up  to  10% 
were  noted,  however.   The  deviations  v/ere  believed  to  come 
from  the  margin  of  error  inherent  in  calculating  the  are^s 
under  the  cross  section  curves. 


TABLE  I 
EFFECTIVE  TARGET  THICKNESS 


GAS 

RUN^ 

e 

Ei(Mev) 

tgff (cm) 

H2 

"2 

«2 
H2 

H2 

Hel 
He2 
He3 
He4 
ne5 

90° 

90° 

120° 

120° 

75° 

55.96 
80.10 
85.90 
63.18 
52.30 

.228±.003 
.234±.003 
.8311.012 
.941±.014 
.6831.010 

"2 
D2 

H2 

D2 

H2 
D2 

"2 
D2 

001 

A  A  T 

004 
004 
006 
006 
009 
009 

90° 

O  A  O 

90° 
90° 
120° 
120° 
120° 
120° 

56.82 

85.82 
70.42 

46.55 

.2741.004 
9  004-  nn  A 

.2251.004 
.2651.004 
.9001.014 
.9441.014 
.9831.015 
1.031.015 

a:  RUNS  001,  004,  006  and  009  were  H2-D2  mixtures  used  by 
Mader  [Ref .4] . 
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Gaby  [Ref.3]  has  used  a  geometrical  approach  to  calculate 
effective  target  thicknesses.   The  geometrical  approach,  how- 
ever, is  very  idealized  and  requires  many  assumptions,  so 
that  it  is  not  as  reliable  as  the  experimental  method. 


B.   THE  SCHWINGER  CORRECTION 

The  Schwinger  radiative  correction  accounts  for  energy 
loss  by  electrons  in  the  field  of  the  target  nucleus,  includ- 
ing nuclear  bremsstrahlung ,  virtual  photon  exchanges,  and 
the  emission  of  photons  below  the  cutoff  energy  of  the 
cross-section  calculations.   The  Schv;inger  correction  is 
much  larger  than  the  Bethe-Heitler  correction  for  the  same 
experiment,  typically  ranging  from  10%  to  20%. 

There  are  several  modifications  to  the  original  correc- 
tion developed  by  Schwinger  in  1949  [Ref.8].   The  version 
used  in  NPGLINAC  calculations  was  put  forth  by  Tsai  in  19  61 
[Ref.9].   The  Tsai  version  of  the  Schwinger  correction  is 


Ks  =  e-«s. 


6g  =  2a 


77 


1/2  In  El  +  1/2  In  E3  -13 
n2AE  AE   12 


in   q2  -1 
m2c2 


137 


where  a  =  fine  structure  constant=   1 
m      =  electron  mass. 


q*-  =  four-momentum  transfer  squared  (in  F  2)  =  q   -  q? . 


and  Ei,E3,  AE  and  n  are  the  same  quantities  that  appeared  in 
the  Bethe-Heitler  correction.  The  quantity  AE  was  chosen  to 
be  at  least  four  half-widths  of  the  peak  of  the  cross  section 
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In  all  the  hydrogen-helium  runs  the  half-v/idth  was  close  to 
0.25  MeV,  so  a  value  of  AE  =  1  Mev  was  used  in  the  calcula- 
tions, not  only  for  the  radiative  corrections  but  for  .. 
deterinining  the  areas  under  the  cross-section  curves  as  well. 
Tsai  also  developed  a  more  complicated  expression  for  the 
Schwinger  correction.    It  is  long  and  difficult  to  evaluate 
except  by  computer.   Mader  has  shown  [Ref.4]  that  for  ener- 
gies obtainable  at  NPGLINAC,  the  more  complicated  Tsai  cor- 
rection did  not  differ  significantly  from  Equation  2-2.   In 
fact,  for  gas  mixture  experiments  v/here  only  the  ratio  of  the 
Schwinger  corrections  for  the  two  gases  was  needed,  the  dif- 
ference  in  the  corrections  calculeited  by  the  com.plicated 
expression  and  Equation  2-2  vras  less  than  0.1%.   Since  most 
of  the  gas  target  experiments  conducted  ar.  NPGLINAC  have  been 
performed  using  gas  mixtures,  it  was  considered  entirely 
adequate  to  use  Equation  2-2  to  calculate  the  Schwinger 
correction  terms. 

III.  HELIUiM  FORM  FACTOR 

A.   INTRODUCTION 

Elastic  electron  scattering  from  the  He^  nucleus  has  been 
performed  several  times,  by  Hofstadter,  et  al .  [Ref s . 10 , 11] , 
by  Erich,  Frank,  Kaas  and  Prange  [Ref. 12],  and  by  Frosch, 
McCarthy,  Rand  and  Yearian  [Ref. 6].   At  low  values  of  q^ ,  the 


•'-  Equation  3-6  in  Ref.  5. 
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charge  form  factor  for  Ke"^  was  found  to  be  in  agreement  with 
the  Gaussian  model: 

F  (q2)  =  exp   [-<r2>q2/6]  (3-1) 

where  <rr'>    is  the  square  of  the  nns  radius  of  the  nuclear 
charge  distribution.   For  He  ,  <r2>l/2  was  determined  to  be 
1.63  F  by  Frank,  Haas  and  Prange  [Ref.l2].   However,  at 
higher  values  of  q2  (q2>6.2  F~2)  a  deviation  from  the  Gaus- 
sian model  v/as  discovered.   In  fact  a  diffraction  zero  was 
found  at  q--  =  10  F~2  [Ref  .6]  .   A  new  expression  for  the  He^ 
form  factor  vras  obtained  on  the  basis  of  the  experimental 
work  of  Frosch,  et  al. 

F(q2)  -   fl  -  (a2q2)n'j  ^^^    (_b2/q2)  (3_2) 

with  n  =  6 

a  =  0. 316  ±  0.001  F 
and  b  =  0.681  ±  0.002  F. 
The  resulting  value  of  <r2>l/2  was  1.68  F, 

For  energies  obtainable  at  NPGLINAC,  it  was  only  possible 
to  study  the  helium  form  factor  at  values  of  q2  less  than 
0.6  F""2^   Therefore  it  was  not  possible  to  distinguish  betv/een 
the  two  models.   A  comparison  V7as  made  to  the  Gaussian  model, 
however,  to  help  determine  the  efficiency  of  NPGLINAC  with 
respect  to  gas  targets  in  general. 
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B .   THEORY 

Linear  accelerators  are  very  useful  in  producing  energe- 
tic electrons  which  can  be  scattered  from  nuclei  in  order  to 
learn  something  about  the  charge  distributions  of  the  differ- 
ent nucleons  and  nuclei.   The  original  theories  of  scattering 
cross  sections  were  developed  by  Rutherford  and  Mott,  and 
assumed  that  nuclei  were  point  particles.   Nuclei  are  not 
point  particles,  however,  and  their  charge  distributions  have 
a  finite  size.   In  the  first  Born  approximation,  the  ratio  of 
an  experimental  cross  section  to  the  Mott  cross  section 
(including  recoil)  is  called  the  charge  form  factor  (F)  and 
is  always  less  than  one  for  values  of  q^>0. 


\  dQ  /  exp/  I,  dil  J 


Mott    =  F^  (3-3) 


where 


(^\  /ge^y   cos  2  2   /  1  \  (3-4) 

\  dQ/  Mott     =\2Eiy  e     \  n  / 


sin  ^  2 


It  can  be  shown  that  the  charge  form  factor  is  just  the 
Fourier  transform  of  the  charge  distribution,  and  is  a 
function  of  the  momentum  transfer: 

F(K)  =    1       (((exp  (-iK-r)  p  (r)d3r 

ge  JJJ  (3-5) 

->  ->-->-> 

where  q  is  related  to  K  by  q  =  fi  K 

For  a  spinless,  spherically  symmetric  nucleus  such  as  He^, 
the  charge  form  factor  can  be  expressed  as 

F(k2)  =1-1  K2<r2>  +  ...   .  (3-6) 

6 
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A  plot  of  F(q  )  vs  g^  should  have  a  y-mtercept  of  1.0  and 

a  slope  of  2^<r2>.   in  this  manner  the  inns  radius  of  the 

6 
charge  distribution  can  be  determined. 


C.   EXPERIMENTAL  METHOD 
1 .   Procedure 

The  experiment  consisted  of  using  NPGLINAC  to  scat- 
ter electrons  from  a  hydrogen-helium  target  at  different 
values  of  q^  and  determiining  the  helium  charge  forro  factor 
from  the  resulting  data.   The  values  of  q^  ranged  from 
0.1  F~2  to  0.55  F~2 .   Appropriate  angles  and  machine  energies 
were  chosen  to  give  the  most  reliable  data,  taking  into  con- 
sideration the  energy  limitations  of  the  machine,  the  length 
of  time  required  for  a  run,  and  the  background  radiation 
levels. 

The  data  for  each  run  consisted  of  a  spectrum  of 
counts  per  MeV  versus  electron  energy  in  MeV.   Spectrometer 
settings  were  such  that  the  elastic  peaks  of  hydrogen  or 
helium  appeared  between  channels  four  and  five  of  the  ten- 
channel  counting  system.   The  total  number  of  points  taken 
during  a  run  was  dependent  on  the  average  energy  separation 
between  counters.   In  addition  to  scanning  the  helium  and 
hydrogen  peaks,  a  number  of  points  were  taken  above  the  helium 
peak  and  between  the  helium  and  hydrogen  peaks  in  order  to 
obtain  measurements  of  the  radiation  background  levels.   The 
integration  of  the  incident  beam  was  chosen  such  that  at 
least  10,000  counts  were  obtained  under  each  oeak.   The  data 
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for  each  point,  including  counts  for  each  of  the  ten  chan- 
nels, integration  capacitance  (yF)  and  voltage,  time  of  inte- 
gration, incident  electron  energy  (MeV)  and  spectrometer 
setting  (MeV)  were  recorded  on  a  teletype  machine. 
2 .   Experimental  Errors 

Because  the  experiment  involved  calculating  the  ratios 
of  areas  under  the  peaks  of  helium  to  hydrogen,  systematic 
errors  in  the  experimental  parameters  cancelled  out.   These 
errors  included  the  deviation  of  the  incident  electron  beam 
energy  from  the  set  value,  which  varied  less  than  0.1%  over 
several  hours;  the  drifting  of  the  magnetic  field  in  the 
spectrometer,  which  was  held  constant  to  within  1  part  in 
10^  by  a  detailed  balancing  circuit;  and  changes  iii  the 
secondary  emission  monitor  (SEM)  efficiency.   While  tne  ab- 
solute SEM  efficiency  was  not  critical  for  evaluating  the 
cross -section  ratios,  it  was  more  important  in  determining 
effective  target  thicknesses  of  individual  gases. 

Statistical  errors  appeared  in  both  the  counting 
rate  and  the  efficiencies  of  the  counters.   The  efficiencies 
of  the  counters  have  been  measured  by  Stewart  [Ref.l3].   With 
a  counting  rate  of  at  least  10,000  counts  per  peak,  the 
statistical  error  was  1.0%  or  less. 

The  largest  error  was  found  to  occur  when  evaluating 
areas  under  the  peaks.   This  was  due  to  difficulties  in 
determining  how  much  of  the  radiative  tail  of  helium  should 
be  subtracted  from  the  hydrogen  peak. 
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D.   DATA  REDUCTION 

1.   Calculation  of  Cross  Sections 

The  raw  data  for  each  run  consisted  of  a  teletype 
printout  of  counts  per  channel;  incident  beam  energy,  spec- 
trometer setting,  and  integration  data.   This  data  was  re- 
duced using  a  computer  program  which  contained  counter 
efficiencies  and  the  energies  "seen"  by  each  of  the  counters. 
The  reduced  data  was  used  to  make  ci  plot  of  energy  versus 
counts  per  microcoulomb .   Two  examples  of  these  plots  are 
shown  in  Figure  2  and  F.igure  3.   They  represent  values  of  q- 
of  0.1031  F~^  and  0.5496  F~^  and  show  the  different  separa- 
tions of  the  peaks  corresponding  to  different  recoil  energies 
of  hydrogen  and  helium. 

The  differential  cross  sections  v;ere  calculated  using 
the  following  relation: 


V  df^  / 


Nsc 


exp    N^N^A^  (3-7) 


where  N    =  number  of  scattered  electrons 

=  Area  under  Peak x  Radiative 

Avg.  Energy  resolution  of  counters   Corrections, 

N-  =  number  of  incident  electrons, 

N^  =  number  of  nuclei  per  unit  area, 
and  A^     =  solid  angle  subtended  by  the  spectrometer 

entrance. 
It  is  also  known  that  the  experimental  cross  section  equals 
the  Mott  cross  section  (cTj^j)  multiplied  by  the  charge  form 
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factor  squared  {Equation  3-3] .   Thus  the  form  factor  for  a 
gas  can  be  obtained  from 

f2  =    N^^ 


NiN^A.Qa^  (3-8) 

In  this  experiment  is  was  desired  to  calculate  the  ratio  of 
the  helium  charge  form  factor  to  the  hydrogen  form  factor. 
Nj^  and  A^  fall  out  of  the  calculation  right  away,  for  they 
are  the  same  for  both  gases. 

F^H      (Nsc'h  (Nt)He  '^M'He 

( 3-9 ) 

The  average  energy  resolution  of  the  counters  is  directly 

proportional  tc  the  energy  of  the  scattering  peak.   In  this 

experiment  also,  the  correction  for  ionization  and  the 

Bethe-Heitler  correction  v/ere  the  same  for  both  gases.   The 

only  radiative  correction  that  appeared  then  was  the  Schv/inger 

correction  K  .   Therefore  the  ratio  of  the  cross  sections 

becomes 

F^He    (A)Hg  (E3)h  (Nt)H  (Om)h  (K^Jh 


f2jj     (A)h   (E3)He  (Nt)He  'OM^He  '^s^He  '^'^°' 

The  areas  under  the  peaks  were  determined  using  a 
three-point  method  of  integration  from  a  prepared  calculator 
program.   The  appropriate  backgrounds  were  then  subtracted 
out.   The  upper  limit  on  the  energy  was  chosen  to  be  where 
the  helium  curve  dropped  to  background  level,  which  usually 


24 


was  about  0.3  MeV  above  the  He  peak.   The  lower  limit  was 
chosen  to  be  about  four  half-widths  below  the  peak  energy, 
which  was  close  to  1.0  MeV  for  all  runs.   The  energy  range 
used  for  the  integrations  was  essentially  the  same  for  both 
the  hydrogen  and  helium  peaks. 

The  gas  concentration  N,  was  calculated  from  the 
known  molar  concentrations  of  H2  and  He.   All  runs  were  con- 
ducted using  a  gas  mixture  that  contained  44.9%  hydrogen  and 
55.1%  helium.   There  are  two  H2  nuclei  for  each  He  nucleus. 
Thus  the  ratio  of  nuclei  concentrations  was 

(N^L  ^  2(44.9)  ^  1.62976 
(Nt)He    (55.1) 

The  Mott  cross  sections  were  computed  from 
Equation  3-4  using  the  desk  calculator,  and  the  Schwinger 
radiative  corrections  were  calculated  from  Equation  2-2.   The 
hydrogen  form  factors  were  obtained  from  a  prepared 
calculator  program.   Finally,  the  helium  charge  form  factors 
were  calculated  from  Equation  3-10. 
2 .   Error  Analysis 

All  errors  in  the  calculation  of  the  helium  charge 
form  factor  were  1.0%  or  less  except  for  the  error  in  eval- 
uating areas  under  the  cross- section  curves.   Specifically, 
it  was  difficult  to  determine  exactly  how  much  to  subtract 
out  from  the  hydrogen  peak  due  to  background  and  the  helium 
radiative  tail.   One  method  used  v/as  to  take  an  average  back- 
ground rate  between  the  two  peaks  and  subtract  it  from  the 
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hydrogen  peak.   At  low  values  of  q^ ,  however,  the  peaks  were 

very  close  together  and  this  made  it  difficult  to  determine 

an  actual  background  rate  [Figure  2].    Another  method  used 

was  to  make  a  least-squares  fit  of  the  helium  tail  from 

halfway  up  the  peak  to  the  start  of  the  hydrogen  peak.   A 

program  was  prepared  to  perform  a  least-squares  fit  of  the 

form  f(x)  =  a  +  bx  +  cx^ , 

where  Xj_  =     1      ,  En  being  the  helium  peak  energy, 
Eq  -  E^ 

and  f  (xj_)  =  counts  per  microcoulomb  corresponding  to  energy 

Ej_.   The  least-squares  fit  worked  well  at  low  values  of  q^ , 

but  at  q^  -  0.55  F~^  the  quadratic  term  caused  the  fit  to 

begin  rising  under  the  hydrogen  peak.   In  this  case  the  fit 

was  reduced  to  f (x)  =  a  +  bx  and  results  were  better. 

The  differences  in  calculations  using  an  average 

background  rate  and  a  least-squares  fit  were  up  to  10%.   The 

method  of  least  squares  resulted  in  form  factor  values  more 

closely  correlated  to  the  expected  results.   Therefore  the 

least-squares  method  was  used  for  all  calculations  in  this 

experiment. 

E.   RESULTS  AND  CONCLUSIONS 

2  ') 

The  experimental  values  of  F(q  )  for  each  value  of  q 

used  are  listed  in  Table  II  and  are  depicted  graphically  in 
Figure  4.   The  error  flags  indicate  the  maximum  statistical 
error  of  1.5%.   The  theoretical  curve  for  a  Gaussian  poten- 
tial is  also  shown  for  comparison.   Although  the  experimental 
values  of  F(q  )  differed  somev/hat  from  the  theoretical  values 
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TABLE  II 
EXPERIMENTAL  He^  CHARGE  FORM  FACTOR 


RUN 

q2(F-2) 

(Fp,^)  exp 

(Fpg)  Gaussian 

1 

.1585 

.982±.015 

.9322 

2 

.3328 

.813±.012 

.8630 

3 

.5496 

.783±.012 

.7840 

4 

.3000 

.815±.012 

.  8756 

5 

.1031   . 

.957±.014 

.9554 
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at  low  values  of  q^ ,  a  least-squares  fit  to  the  form 
y  =  ae^^  produced  a  curve  that  nearly  coincided  with  the 
theoretical  curve.   The  resulting  r.ias  radius  of  the  He"^ 
nucleus,  obtained  from  the  least-squares  fit,  was  <r2>l/2 
=  1.78  F,  larger  than  the  expected  value  by  9.2%.   The  dis- 
crepancies were  believed  to  be  due  mainly  to  the  difficul- 
ties encountered  when  evaluating  areas  under  the  hydrogen 
peaks.   At  low  values  of  q^  the  separation  between  peaks  was 
not  wide  enough  to  obtain  an  accurate  determination  of  back- 
ground and  helium  radiative  tails.   It  was  noted  also  that 
for  q2  =  .1585  F~2  the  background  level  increased  signifi- 
cantly from  the  helium  peak  to  the  hydrogen  peak,  causing  the 
experim.ental  form  factor  to  be  larger  than  predicted.   At 
q^'  =  .55  however,  tliere  was  good  agreement  betv/een  experiment 
and  theory.   This  was  believed  due  to  adequate  separation 
between  peaks  for  an  accurate  background  determination,  and 
to  a  fairly  constant  background  rate  throughout  the  energy 
range  of  the  calculations. 

It  is  believed  that  better  agreement  between  helium 
form-factor  experiments  at  NPGLINAC  and  those  at  other  lab- 
oratories can  be  obtained  by  using  pure  gas  targets  instead 
of  mixtures.   This  would  ease  the  burden  of  attempting  to 
calculate  the  proper  background  rates. 
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APPENDIX  A 
FREQUENTLY  USED  CONSTANTS 

The  following  physical  constants  appeared  several  times 
in  cross-section  calculations  and  are  listed  here  for  future 
reference: 

e^  =  0.144  X  10"-^^  MeV  -  cm 

(Mc2)^  =  938.68  MeV 

(Mc2)j^^  =  3728.01  MeV 


m^c    386.12  F 


2a  =  _ 

TT     0.46455  X  10  "^ 


]  CTTT^  = 


,5  ,^      2 


erg     0.68947  x  10^  in. 

lb 


RT  =  0.643041  x  10^^  erg    at  T  =  77.32°  K 

mole 

1  dyne  =  0.14504  x  10"^   lb 


,2 


cm^  in. 2 
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APPENDIX  B 
VlRIAL  COEFFICIENTS 

The  theoretical  development  of  the  virial  theorem  and 
the  virial  equation  of  state  is  described  in  Ref.l4  and 
Ref.l5.   The  virial  equation  of  state  is  an  alteration  to 
the  ideal  gas  law  (pV  =  nRT)  and  more  accurately  describes 
the  properties  of  real  gases.   The  virial  equation  of  state 
appears  in  two  forms,  a  power  series  expansion  in  the  molar 
volume  and  a  power  series  expansion  in  the  pressure.   Since 
pressure  is  a  directly  measureable  quantity  for  gas  target 
runs  at  NPGLINAC,  the  latter  form  is  more  convenient  to  use: 


£V  =-  1  +  b'  (T)p  +  C"(T)p^  +  ...    .  (B-1) 

RT 


V   =  V   ,  the  molar  volume,  where  Nq  is  Avogadro ' s  number 
No 

B"(T)  -  B/RT  and  C" (T)  =  (C-B^)  /  (RT)^, 


where  B  and  C  are  the  second  and  third  virial  coefficients 
and  are  functions  of  temperature  only.   For  the  gases  current- 
ly being  studied  at  Monterey  (hydrogen,  deuterium  and  helium) 
the  virial  coefficients  were  calculated  from  the  Lennard-Jones 
potential  for  non-polar  molecules.   Preliminary  calculations 
showed  that  for  a  pure  gas  target,  the  values  of  B''(T)  were 

of  the  order  of  10    in.  /lb  and  values  of  C  (T)  were  of  the 

-7     2 
order  of  10    in.  /lb.   The  terms  involving  higher  orders  of 

p  were  dropped  with  no  significant  changes  in  results. 
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Evaluation  of  B(T)  and  C(T)  for  pure  gases  is  straight- 
forward, using  the  relations 


B   =    b^    B*(T*)     and    C    =    b'^    C*(T*). 


2  *? 

b^   =    2/3TrN.    a      cm-^/mole,    and   T*    =   kT/e , 


'0  ~  ^/-^"^^o 

the  reduced  temperature,  o,    c/k  and  bQ  are  force  constants 
for  the  Lennard-Jones  potential  and  are  tabulated  for  sever- 
al gases  [Ref .15] ,  including  those  in  use  at  Monterey.   Two 
values  for  each  parameter  are  quoted,  one  determined  classi- 
cally and  the  other  quantum  mechanically.   The  quantum- 
mechanical  parameters  are  more  accurate  in  principle  and  were 
used  in  all  calculations. 

Once  T*  has  been  determined,  B*  and  C*  can  be  extracted 
from  a  table  [Ref. 15].   For  values  of  T*  not  exactly  equal 
to  the  tabulated  values  of  B*  and  C* ,  a  linear  interpolation 
is  necessary  to  obtain  correct  values  of  B*  and  C*.   Proper 
substitution  yields  the  quantities  B^(T)  and  C^  (T)  ,  vrhich  can 
then  be  entered  into  the  equation  of  state.   All  parameters 
are  then  known  except  V,  v/hich  is  related  to  the  desired 

parameter  p (density)  by  V  =  W,  where  W  is  the  molecular 

P 

weight.   Solving  for  the  density  yields  the  following 
formula: 


p   =  I  =  "p 


^    RT  [  1  +  B"(T)p  +  C''(T)p2] 

(B-2) 
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The  calculations  of  the  virial  coefficients  Bj^j_^  and 
Cj^j_-^    for  a  gas  mixture  are  slightly  more  complicated.  In  fact, 
^ix  ^^^  ^°^  ^^  calculated  exactly  so  it  is  dropped  from  the 
expansion.   In  general,  ^mix  ^^^  ^  5^^  mixture  is  given  by 
n    n 

^mix  ^      /       /   ^i  ^j  ^ij  '  (B-3) 


i=l  j=l 

where  n  =  number  of  constituent  gases  in  the  mixture  and 
Xj_  =  molar  fraction  of  the  ith  gas  in  the  mixture. 
B- •  (i  7^  j)  is  a  theoretical  virial  coefficient  for  the  gas 
mixture,  v/here  T-  •*  and  bg  are  determined  from  the  nev;  para- 


meters e^-i  =  /e'.    £~  and  o ^  j,    =    1/2  (a^  +    o-)  ,    obtained  from 
the  force  constants  of  the  ith  and  jth  gases  in  the  mixture. 
B- •*  is  extracted  from  the  table  in  the  same  manner  as  be- 
fore, only  the  table  is  entered  with  the  value 
•^ii*  ~  kT/ej_-j  .   B-;j  is  merely  the  second  virial  coefficient 
for  the  pure  jth  gas  in  the  mixture. 

For  typically  used  gas  mixtures  (n  =  2) ,  the  second 
virial  coefficient  becomes 

Bmix  =  ^l    Bll  +  2X1  X2  B12  +  X^  B22.  (B-4) 

The  temperature  of  all  gases  and  gas  mixtures  used  at 
NPGLINAC  was  the  temperature  of  the  liquid  nitrogen  that  was 
used  as  a  coolant  (77.32°  K) .   Variations  from  this  tempera- 
ture were  small  during  all  runs,  and  it  was  assumed  constant 
for  all  calculations.   Figure  5  depicts  the  temperature 
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bridge  control  mechanism.   Values  of  the  virial  coefficients 
for  the  gases  used  are  listed  in  Table  III.   Table  IV  shov/s 
values  of  gas  densities  for  typical  pressures  calculated 
from  the  virial  equation  of  state. 
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TABLE  III 
VIRIAL  COEFFICIENTS  FOR  SEVERAL  GASES  AT  T  =  77.32°  K 


> 

GAS 

B"(T)  (in.2/ib) 

C"(T) (in.^/lb^) 

«2 

-.18966  X  10"^ 

+  .12928  X  10"'^ 

°2 

-.18966  X  10"^ 

+  .12928  X  10""^ 

48  8%  H'^- 

« 

-.18944  X  10"^ 

51.2%  D2 

mixture 

\ 
He 

+.89735  X  10"^ 

+  .71852  X  10""^ 

44.9%  H2- 
55.1%  He 

+.43683  X  10~5 

mixture 
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TABLE  IV 

DENSITIES  OF  SEVERAL  GASES  AT  CONSTANT 
TEMPERATURE  (77.32°  K)  FOR  VARIOUS  PRESSURES 


Gas 

p(lb/in.2) 

p  (Qin/cm^)    X 

10-2 

H^  (in 

140 

.310906 

148 

.329184 

149 

.331473 

H2-D2 

150 

.333699 

151 

.336053 

mixture) 

152 

.338344 

160 

.356709 

140 

.621334 

D9  (in 

148 

.657862 

149 

.662437 

H2-D2 

150 

.666885 

151 

.671591 

mixture) 

152 

.676170 

160 

.712871 

140 

.302475 

H2  (in 

148 

.319748 

149 

.321907 

Ho-He 

150 

.324066 

^ 

151 

.326225 

mixture) 

152 

.328384 

160 

.345655 

140 

.600643 

He  (in 

148 

.634943 

149 

.639231 

H2-He 

150 

.643518 

151 

.647805 

mixture) 

152 

.652093 

160 

.686389 
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TABLE  IV  (Cont.) 

DENSITIES  OF  SEVERAL  GASES  AT  CONSTANT 
TEMPERATURE  (77.32°  K)  FOR  VARIOUS  PRESSURES 


Gas 

p(lb/in.^) 

p  (gm/cm-^)  x 

10~2 

140 

.310834 

H2 

148 

.329099 

149. 

.331386 

(pure) 

150 

.333839 

151 

.335963 

152 

.338252 

160 

.356600 

140 

.621204 

D2 

148 

.657694 

.  149 

.662264 

(pure) 

150 

.667166 

151 

.671411 

152 

.675986 

160 

.712654 

140 

.593471 

He 

148 

.626929 

149 

.631108 

(pure) 

150 

.635286 

151 

.639463 

152 

.643640 

160 

.677024 
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APPENDIX  C 
GAS  TARGET  EFFECTIVE  THICKNESS 

The  thickness  of  a  solid  target  is  directly  measureable 
and  poses  no  problem  in  a  cross-section  calculation.   How- 
ever, a  gas  target  is  different.   The  pressurized  gas  is  con- 
tained in  a  cylindrical  chamber  two  or  three  inches  in  . 
diameter,  but  only  those  electrons  scattered  from  a  small 
volume  of  the  gas  make  their  way  through  the  spectrometer 
entrance  slit.   The  calculation  of  the  effective  thickness 
of  a  gas  target  is  detailed,  and  is  reviewed  here  to  show 
what  was  done  for  this  experiment,  and  to  clarify  the 
procedure  for  future  reference. 

The  experimental  cross  section  is  given  by  Equation  3-7. 
The  target  thickness  t  is  located  in  the  term  N^ ,  the  number 
of  target  nuclei  per  unit  area. 


Nt  =    PtNQ 


(C-1) 


where  A  is  the  atomic  weight,  p  is  the  gas  density  determined 
using  the  virial  equation  of  state,  and  Nq  is  Avogadro's 
number.   Substituting  Equation  C-1  into  Equation  3-7  and 
solving  for  t  gives  the  follov^ing  expression: 


t  = 


N 


sc 


N.  A^ 
exp   1 


A 


cN, 


(C-2) 
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N    is  the  number  of  scattered  electrons  and  is  equal  to 
the  integrated  area  (counts-MeV)  divided  by  the  average  reso- 
lution of  the  counters  AE^  (MeV). .   AE^.  can  be  calculated  from 
the  dispersion  formula, 

AE-p  _  Ap  _  a/r 


E     p      D  (C-3) 

where  a  =  3/16" 

r  =  16" 

and  D  =  Dispersion  factor  =  3.92. 

Thus  E^=  0.00299  E. 

The  solid  angle  A^  =  Area  of  spectrometer  entrance  slit/r^ 

=  (5/16) (1  1/2) 
162 

=  0.183  X  10~2    sr.. 

The  number  of  incident  electrons  is  determined  from  the  total 
charge  accumulated  in  the  SEM  and  the  SEM  efficiency: 

Ni  =    Q 

e  EsEM  (C-4^ 

e  is  the  charge  of  one  electron.   The  SEM  efficiency  v/as 

assumed  constant  at  6.1%. 

Now  the  actual  integrations  calculated  from  the  reduced 

data  were  in  units  of  counts  per  yC  -  MeV.   This  causes  the 

Q  in  N-  to  drop  out  when  N   /Nj  is  calculated: 
1        ■^  sc   1 

Nsc  =  ^^e^-e-Eggj^ 


N^        AE^  (C-5) 
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The  experimental  cross  section  is  evaluated  using  Equation 
3-3.   Thus  the  final  form  of  the  expression  for  effective 
target  thickness  becomes: 


t  =  Area^e'Eg^j.^'A 
aj^-AEj--AQ«p'NQ 

=  (0.29605  X  10-32)  f    AreavA)     cm.  (C-6) 
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